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Although the Kondo effect and the Kondo ground state of a magnetic impurity have 
been investigated for more than forty years it was until recently difficult if not impos- 
i sible to calculate spatial properties of the ground state. In particular the calculation 

of the spatial distribution of the so-called Kondo cloud or even its existence have been 
elusive. In recent years a new method has been introduced to investigate the properties 
(~~~. | of magnetic impurities, the FAIR method, where the abbreviation stands for Friedel 

Artificially Inserted Resonance. The FAIR solution of the Friedel-Anderson and the 
. Kondo impurity problems consists of only eight or four Slater states. Because of its 

compactness the spatial electron density and polarization can be easily calculated. In 
■ this article a short review of the method is given. A comparison with results from 

0^ . the large N-approximation, the Numerical Renormalization Group theory and other 

methods shows excellent agreement. The FAIR solution yields (for the first time) the 
# ^ ' electronic polarization in the Kondo cloud. 
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1 Introduction 



The properties of magnetic impurities in a metal is one of the most intensively studied 
problems in solid state physics. The work of Friedel [1] and Anderson [2] laid the foun- 
dation to understand why some transition-metal impurities form a local magnetic moment 
while others don't. Kondo [3] showed that multiple scattering of conduction electrons by a 
magnetic impurity yields a divergent contribution to the resistance in perturbation theory. 
Kondo's paper stimulated a large body of theoretical and experimental work which changed 
our understanding of d- and f-impurities completely (see for example [4], [5], [6], [7], [8], [9], 
[10], [11], [12], [13]). A large number of sophisticated methods were applied in the following 
three decades to better understand and solve the Kondo and Friedel-Anderson problems. 
In particular, it was shown that at zero temperature the Friedel-Anderson impurity is in 
a non- magnetic state. To name a few of these methods: scaling [14], renormalization [15], 
[16], [17], [18] Fermi-liquid theory [19], [20], slave-bosons (see for example [21]), large-spin 
limit [22], [23]. After decades of research exact solutions of the Kondo and Friedel-Anderson 
impurities were derived with help of the Bethe-ansatz [24], [25], [26], representing a magnif- 
icent theoretical achievement. The experimental and theoretical progress has been collected 
in a large number of review articles [7], [8], [9], [10], [11], [12], [13], [15], [20], [21], [23], [24], 
[25], [26], [27]. 

The exact theory of the Bethe ansatz is such a complex theory that only a limited number 
of parameters can be calculated. For the majority of practical problems one uses the numer- 
ical renormalization group (NRG) theory and the large-spin (large N) method. Recently 
the author introduced another approximate solution for the Friedel-Anderson (FA) [28], [29] 
and the Kondo impurity [30] , the FAIR method. The FAIR solution consists of only four to 
eight Slater states and is therefore very compact. It is well suited to calculate in particular 
spatial properties of the Kondo ground state. It yields the first quantitative calculations 
of the Kondo cloud [31]. There are very few spatial properties of the Kondo ground state 
calculated with other theoretical approaches. One example is the NRG calculation for the 
Friedel oscillations in the vicinity of a Kondo impurity [32]. A calculation of the Friedel 
oscillations with the FAIR method yields good agreement with the NRG results [33] . In this 
short review the FAIR method will be introduced and some of the results presented. The 
FAIR method uses Wilson states [15] which replace a complete conduction electron band by 
a relatively small number of states which carry the full interaction with the impurity. The 
Wilson states are sketched in appendix A. 

2 II The FAIR Method 

2.1 The artificial Friedel resonance state 

We consider the Hamiltonian of a band with a finite number N of non-interacting electron 

states 

#0 = Y£=1 £ v<t,Cv 
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The <4 are the creation operators of the band. In the following the states such as <4$o are 
represented and addressed by their creation operators cj, suppressing the vacuum states $o- 
From these band states a new (arbitrary) state aj is composed 



N 



a o = $^«o c t (!) 



u=l 



In the next step an intermediate basis | a l| can be constructed numerically where the 

additional (N — 1) states a\ are orthonormal to each other and to a\. In this basis the 
Hamiltonian H is given by an N x N matrix with the elements (Ho)^ where (H ) 00 is at the 
left upper corner. In the final step the (N — 1) x (N — 1) sub-matrix of (Hq)^ for i,j>0 

is diagonalized. This yields the new basis | a l| — { a o> a L a L ■■■> a jv-ij which is uniquely 
determined by the state a\. In this basis the s-band Hamiltonian has the form 



N N 



H = E i a \ a i + E ala + ^ V fr (i) ^aja, + a\a^j (2) 

u=l v=\ 

One recognizes that this Hamiltonian is analogous to a Friedel Hamiltonian where a\ repre- 
sents an artificial Friedel resonance. Therefore this state is called a FAIR state for Friedel 
Artificially Inserted Resonance state. 

It has to be emphasized that the FAIR state a\ can have any composition of the basis 
states cj,. Therefore it can be adjusted to a given problem without any restriction. This 
gives the FAIR method its adaptability. 



2.2 The Friedel resonance 

As an example let us consider the simple Friedel resonance Hamiltonian Hp T . 

N N 

H FR = J2 £ » c " c » + E ^ d + J2 Vsd ^ + ( 3 ) 

Since Hfr does not depend on the spin the latter will be ignored. The first term is the 
conduction band Hamiltonian Hq, the second term gives the energy of the d (resonance) 
state of the impurity with being its creation operator. The last term represents the 
interaction between the d state and the conduction electrons. 

There exists a FAIR state aj and a FAIR basis | a j| so that the n-electron ground state 
of the Friedel Hamiltonian is exactly given by 

y Fr =(Aal + BS Ul«I^o (4) 

v ' i=i 
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Here A and B are coefficients which fulfill the condition A 2 + B 2 = 1. Actually this exact 
form of the ground state of the Friedel impurity can be understood without any analytic or 
numerical calculation [34]. This is shown in the appendix D. In ref. [35] it was discovered 
by a variation of a\ minimizing the ground state energy of the state (4) with respect to the 
Friedel Hamiltonian (3). The state a\ determines all the other basis states a\ uniquely. Since 

the new basis nas the same number of states as the original basis {cj,} the construction 

of the basis is onr y possible if the number N of basis states is small. For N ~ 10 23 it 

would be hard to construct the orthonormal sub-diagonal basis j°! j • Wilson has shown in 

his Kondo paper [15] how to construct a finite basis {cj,} which preserves the full interaction 
with the impurity. The Wilson states are discussed in the appendix A. 

The example of the Friedel Hamiltonian shows the simplicity and effectiveness of the 
FAIR method. It can be applied to treat the Friedel- Anderson and the Kondo impurity. 



3 The Friedel- Anderson impurity 

The FA-Hamiltonian consists of the Friedel Hamiltonian (3) for both spins plus a Coulomb 
term of the form H c = Un^n^. 

{N N ~\ 

^2 e v c\ c c va + E d d\d a + ^ V sd (v) [d\c va + S vc d c \ \ + U n d1 n dl (5) 
v=\ v=\ ) 

3.1 The magnetic state 

In the early years (before the Kondo paper) it was the goal to calculate (and measure) 
the magnetic moment of a d- or f-impurity. After the discovery of the Kondo effect and 
after Schrieffer and Wolff [36] transformed the FA-Hamiltonian into a Kondo Hamiltonian it 
became clear that the ground state of the FA impurity is non-magnetic. Then the calculation 
of the magnetic moment was often considered as irrelevant, even heresy. The paper by 
Krishna- murthy, Wilkins, and Wilson [17] clarified the role of the local magnetic moment in 
the FA-impurity. KWW performed a numerical renormalization a la Wilson [15] for the FA- 
Hamiltonian. They demonstrated that the renomalization-group flow diagram showed very 
different flows from the free-orbital fixed point Hp Q to the strong coupling fixed point Hg C 
(see Fig.l). For sufficiently large Coulomb repulsion (when U » T = np \V sd \ 2 ) the flow of 
their Hamiltonian Hn passed close to the (unstable) fixed point H£ M for a local moment. 
This means that under these conditions the impurity assumed first a magnetic moment 
when the temperature is lowered. After passing the fixed point for the local moment Hl M 
the renormalization flow is essentially the same as for a Kondo Hamiltonian (where a local 
moment is the starting point). 



4 



A 




Fig. 1: Schematic 
renormalizat ion-group flow 
diagram after ref.[17]. 



With decreasing ratio of U /V the flow path passes less and less close to Hl M . This means 
that the size of the local moment decreases until there is no longer a local moment formed. 
(The flow of the susceptibility indicates this behavior). 

At the end point of the renormalization the system approaches the strong coupling fixed 
point Hg C and shows the universal behavior of the Kondo ground state. Nevertheless the 
ground-state wave functions are quite different for small and large ratios of U/F because the 
size of the magnetic moment is engraved into the wave function. 

Let us first consider the local moment state of the FA-impurity. This state is a ground 
state if one applies a magnetic field which is strong enough to suppress the Kondo ground 
state. Within the FAIR approach the (potentially) magnetic solution has the form 



n—\ n—1 

^ms = U^aJA + Kd<iA + a^4 6 !h + 4w44l II 4 T II ( 6 ) 



i=l i=l 



where | a l| an d j^lj are two (different) FAIR bases of the iV-dimensional Hilbert space. 

The state (6) opens a wide playing field for optimizing the solution: (i) The FAIR states aL 

and can be individually optimized, each one defining a whole basis | a l| > j^l} anc ^ ^ e 
coefficients A a ^, A a ^, A^^, A^d can be optimized fulfilling only the normalization condition 
A\ b + A a d + A\ b + A\ d = 1. Since the relative size of the coefficients A a ^, A a ^, Ad y b and A^d 
is not restricted this solution describes correlation effects well. The optimization procedure 
is described in detail in the appendix B. 

Fig. 2 shows the structure of the four Slater states of $ms graphically. The FAIR states 
al and b\ are imbedded in the spin-up and down bands while the 4 an d d\ states are shown 
on the left and right side of the bands. Full circles represent occupied states. 
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Fig. 2: The composition of the magnetic state ^ms is shown. It consists 
of four Slater states. Each Slater state has a half-filled spin-up and 
down band, two FAIR states (circles in within the bands) and two d-states 
(circles on the left and right of the band). Full black circles represent 
occupied states and white circles represent empty states. The band at 
the right with the half-filled circles symbolizes the magnetic solution with 
four Slater states. 

Fig. 3a shows the magnetic moment as a function of U for the mean-field solution and the 
magnetic state ^ms- I n Fig. 3b the ground-state energies of the mean-field solution and the 
magnetic state are compared. The magnetic FAIR solution has a considerably lower energy 
expectation value than the mean-field solution. More importantly it increases the critical 
value of U for the formation of a magnetic moment by almost a factor two (compared with 
the mean-field solution). The mean-field approximation is still used in combination with 
spin-density functional theory (SDFT) to calculate the magnetic moment of impurities [37], 
[38], [39], [40], [41]. A combination between SDFT and the FAIR solution appears to be 
very desirable. 



0.2 



Mean 
fieli 



N=40 
E a =-U oo|i| /2 

|V /=0.1 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 



-0.10 



-0.15 



-0.25 



Mean field 



Magnetic state 



Singlet state 



N=40 
|V/=0.1 
E =-U 72 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Fig. 3a: The magnetic moment as a function of the Coulomb exchange energy U, 
using the mean-field solution and the magnetic FAIR solution *&ms- 
Fig.3b: The ground-state energies of the mean-field solution, the magnetic and 
the singlet FAIR solution 
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3.2 The singlet state 

In hindsight it is quite natural that the magnetic state with its broken symmetry is not the 
ground state. By reversing all spins one obtains a new state with the same energy. In Fig. 4 
the energy of the magnetic state is plotted as a function of the magnetic moment. 




Fig.4: The energy of the FA impurity as a function 
of the magnetic moment. 

This is a situation similar to an atom in a double well potential. In the ground state the 
atom is in a symmetric superposition of the wave functions in the two wells. In analogy one 
can construct the singlet ground state ty$s of the FA-Hamiltonian. This state is obtained 
by reversing all spins in (6) and combining the two wave functions 



^ss = (T) + ^ms (I) 



+ 



n— 1 n— 1 



i=i i=i 

n— 1 n— 1 



Aa^aly + A a4 d\a} ol + A^&^ciJ + A d4 d\d\ [J 6^ JJ aj^o 



i=l i=l 



(7) 



Again one has to optimize a , £»q and all the coefficients. It is remarkable that the 



composition of the FAIR states changes dramatically for small energies. The ground-state 
energy of the singlet FAIR solution lies considerably below that of the magnetic FAIR solu- 
tion ^ms (see Fig. 3b). One can compare this ground-state energy with a set of numerical 
calculations by Gunnarsson and Schoenhammer [42]. They applied the large Nf method 
to the (non-degenerate) FA-Hamiltonian (spin 1/2) for a finite Coulomb interaction and 
included double occupancy of the impurity level. They calculated the ground-state energy 
in the l/iVy-expansion up to the order (1/Nf) 2 which includes more than 10 7 basis states. 
For the s-d-hopping transition they used an elliptic form. With the following parameters: 
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band width B GS = 6eV, Coulomb energy Ugs = 5eV, d-state energy E d ^s = —2.5eV they 
performed two calculations, one for s-d coupling Vqs = leV and another for Vqs = 2eV. The 
table compares the ground-state energies and the occupation for of the d-states {do, d±, d 2 ) 
obtained by GS and the FAIR for Vqs — leV^ and 2eV. Not only the ground-state energies 
but also the occupation of the d-states agree remarkably well. 



V GS = leV 



V GS = 2eV 



states 


E [eV] 


do 


di 


d 2 


no. of coeff. 


GS 


-0.245 


0.034 


0.931 


0.034 


>10 7 


FAIR 


-0.239 


0.035 


0.931 


0.034 


80 




states 


E [eV] 


do 


d\ 


d 2 


no. of coeff. 


GS 


-1.217 


0.137 


0.732 


0.132 


> 10 7 


FAIR 


-1.234 


0.140 


0.722 


0.138 


80 



Table Ia,b: The ground-state energy E and the occupations do, d\, d 2 
of the d-states with 0,1 or 2 electrons. 

It is worthwhile to remember that the FAIR solution is completely determined by the 
two FAIR states a\ and b\, i.e. by 2 x 40 amplitudes for a typical value of iV = 40. On the 
other hand the large N calculation describes the ground state by more than 10 7 parameter, 
i.e., amplitudes of Slater states. 



4 The Kondo Impurity 

The Kondo Hamiltonian is a limiting case of the FA-Hamiltonian [36]. It applies when the 
exchange energy U approaches infinity while the d-state energy approaches — oo, for example 
E d = —U/2. Then the d-state is always singly occupied, either with spin up or down. The 
interaction between the spin s of a conduction electron and the spin S of the impurity can 
be expressed in the form 2Js • S where J > 0. 

In this case the ansatz for the compact FAIR-solution can be obtained from equ. (7). 

The coefficients A a h , A a ^ b ; A d d , A d d have to vanish because there is only single occupancy of 
the d-state in the Kondo solution. This yields 



n— 1 n— 1 



K 



AaAA + A ^ d \ b ii] n 4 n b i®° ( § ) 



i=i i=i 

n—l n—1 



+ 



Aa^a^ + A djb b^d\ J] h n II °a $0 



1=1 1=1 



Again one can optimize the localized states a\ and b\ and the coefficients A a;d , A a ^ d , 



A d< b, A dj b- If one arranges the spin up to the left and spin down to the right in all components 
of (7) then one obtains in the ground state A a ^ d = A a>d and A d ^ = A d ^. Our group calculated 
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the total spin of this state (for J = 0.1) and obtained for the expectation value of (S 2 ) = 
(( s d + J2i s i) 2 ) the value 0.04 in the ground state [30]. For the first excited state one obtains 
(S 2 ) = 1.99. This means that the ground state is essentially a singlet state ((S 2 ) = 0) and 
the first excited state a triplet state ((S 2 ) = 2). 

In the Kondo effect one is generally not so much interested in the ground-state energy but 
in the so-called Kondo energy. This is, for example, the energy difference between the triplet 
and singlet states. The logarithm of this excitation energy is plotted in Fig. 5 as a function 
of 1/ (2Jp ) as the full circles (po is the density of states). The straight line corresponds 
to AE ~ 5Dcxp [—1/ (2Jp )]. This is the unrelaxed singlet-triplet energy which uses for 

the triplet state the same bases | a j| an d j^lj as m the singlet state. One can derive a 
relaxed triplet state by the following trick. In the triplet state the coefficients are related by 
A a4 = -A atd and A djb = -A djb . If one replaces A ajd , A d)b from the start by -A a4 , -A djb and 
optimizes the energy then one obtains the relaxed triplet energy. The difference between 
this energy and the singlet energy yields the relaxed excitation energy AE st . (Since these 
are two independent calculations they have to be performed with an absolute accuracy of 
10~ 10 ). This relaxed excitation energy is plotted in Fig.5 as stars. The stars lie between 
two theoretical curves: (i) AE st = Dexp [—1/ (2Jp )] , given by the dashed curve and (ii) 
AE st = y / 2Jpo_D exp [—1/ (2Jp )], given by the dotted curve. Both expressions are given in 
the literature as approximate values for the Kondo temperature ksTx- The numerical values 
lie closer to the second expression. Therefore the relaxed singlet-triplet excitation energy 
corresponds closely to the Kondo energy and confirms that the FAIR method represents the 
physics of the Kondo impurity accurately. 
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Fig. 5: The energy difference between the singlet and triplet 
states. The full circles represent the unrelaxed singlet-triplet 
excitation energy AE st while the stars yield the relaxed singlet-triplet 
excitation energy AE* t (see text). The dashed and dotted curves are 
theoretical expressions for the Kondo energy. 

The FAIR method yields both energies, the ground-state and the singlet-triplet excita- 
tion energy, with good accuracy although the two energies can differ by a factor of thousand. 

5 Real Space Properties 

Since the compositions of the magnetic state and the singlet state are explicitly known and 
consist only of a few Slater states it is straight forward to calculated the electron density 
and spin polarization of the different states. The details of the calculation are described in 
ref. [31]. 

In Wilson's approach the wave number k is given in units of the Fermi wave number fop. 
Therefore it is convenient to measure real space distances £ in units of i.e, half the 

Fermi wave length. (In this unit the wave length of the Friedel and RKKY oscillations is 

"1"). 
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The density of the Wilson states ip v (£) is given by 



(for v < N/2). The main contribution to the density p° of the state ip v lies roughly in the 
region |£| < 2 U+2 (in units of Xp/2). The different ip v (£) have very different electron densities 
and vary by roughly a factor of 2 N I 2 (which is generally larger than 10 6 ). Therefore it is 
useful to calculate the integrated electron density on a logarithmic scale. 

First we discuss the magnetic state whose wave function is given in equ.(6). 

5.1 The Magnetic State 

The magnetic state ^ms is the building block of the singlet state. Its multi-electron state 
is built from four Slater states and shown in equ. (6). The electron system has already a 
finite density without the d-impurity. Therefore it is useful to calculate the change of the 
(integrated) densities for spin-up and down conduction electrons due to the d-impurity. 

In Fig. 5 these (integrated) densities as well as their sum and difference (total density 
and polarization) are plotted for the parameters E d = —0.5, \V® d \ 2 = 0.04, U = 1 and 
N = 50. With these parameters the impurity has a well developed magnetic moment of 
\i = 0.93/zb- The occupation of the different components is A ab = 0.0294, A ad = 0.0057, 
A db = 0.9355 and A dd = 0.0294. This means that 93.6% of the densities is due to the Slater 

state ^d,b = d\boi YYi=i a \-\ Y\dZ\ The abscissa is the logarithm (using the basis 2) of 

f = 2x/\ F . 

The region beyond £ = 2 20 corresponds to the rim or surface of the sample and is discussed 
below. One recognizes that there is only a negligible polarization of the electron gas in the 
vicinity of the impurity. The important result of Fig. 5 is that there is no polarization cloud 
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around the magnetic state of the impurity. 
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Fig.5: The net integrated density f Q r p(£) <i£ of the s-electron within a distance r 
from the impurity for spin up and down, as well as total density and spin 
polarization. The magnetic moment of the impurity is 0.93^. 

If one looks at the magnetic state, in particular the dominant component ^>d,b = d^b^ YYi=i a \-\ YYi=i 
then one realizes that the total magnetic moment of all the conduction electrons is equal 
to —lfiB (cancelling the moment of +l/iB of the d-electron which is not shown in Fig.5). 
So how can the state $?ms have a finite magnetic moment. The answer is given by Fig.5. 
The moment /ie of the s-electrons is pushed towards the surface of the sample which is at 
the largest radius used for the Wilson states, i.e. 2 N I 2 . This explains the change of the 
integrated densities and polarization at I ~ 22 from zero to —1. If one increases the number 
iV of Wilson states by AiV then the transition is shifted by Al = AN/2. 



5.2 The Kondo cloud 

One of the most controversial aspects of the Kondo ground state is the so-called Kondo cloud 
within the radius tk where rx is called the Kondo length 

}%vf de/dk , . 

r K = = (9) 

(ek= Kondo energy, vp = Fermi velocity of the s-electrons). For a linear dispersion relation 
this yields in reduced units = 1/ (ttek) where in this relation the Kondo energy is given 
in units of the half-band width. 
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The idea is to divide the ground state ^ k of a Kondo impurity into two parts with 
opposite d-spins. The proponents of the Kondo cloud argue that in each component there 
is an s-electon cloud within the Kondo sphere which compensates the d-spin. An important 
assumption of the Kondo-cloud proponents is that, above the Kondo temperature, the bond 
is broken and this screening cloud evaporates from the Kondo sphere. 

In the 1970's Slichter and co-workers [43] investigated Cu samples with dilute Fe-Kondo 
impurities by means of nuclear magnetic resonance. They did not detect any Kondo cloud. 
In a number of recent theoretical papers [44], [45], [46], [47] the argument is made that the 
old NMR experiments could not possibly have detected the screening electron because of the 
large volume of the Kondo sphere yielding a polarization of less than 10 -8 electron spins per 
host atom. 

The FAIR solution of the singlet ground state is well suited to determine the electron 
density and polarization in real space [31]. To the knowledge of the author this is the first 
detailed calculation of the Kondo cloud. 

In the following analysis of the singlet state ^ ss the same parameters are used as for 
the magnetic state in Fig. 5: E d = —0.5, |V,° | 2 = 0.04, U = 1. This yields the following 
squared amplitudes: A^~ b = 0.0146, A^~ d = 0.0028, A^ b = 0.4629 and ~A\ = 0.0146. (The 

A x>y amplitudes are identical). These occupations are very close to half the values of the 
magnetic state {A\ s = 0.0294, A\ d = 0.0057, A\ s = 0.9355 and A\ d = 0.0294). This 



means that ^ S s is given in good approximation as ^ S s ~ (l/\/2~) ^ms (T) + ^ms(1) 
The two magnetic states with d-spin up and down are robust building blocks of the singlet 
state. (However, there are subtle changes in the FAIR states which will be discussed below). 
Therefore the spin polarization of one of the magnetic components, for example of ^ms (t)> 
would be of interest. 

In Fig. 6 the integrated densities of spin up and down electrons, their sum and difference 
(the polarization) are plotted versus the distance from the magnetic impurity (on a logarith- 
mic scale). One recognizes that now one has considerable contributions to the integrated 
net densities of both spins. The polarization of the two contributions is no longer zero 
but reaches a value of —0.46 at a distance of r « 2 11 " 6 . Since the magnetic state $ms (T) 
with net d-spin up has only a weight of about 1/2 it contributes an effective <ij-moment of 
0.93/2 w 0.46. Therefore this d-spin is well compensated by the polarization of the s-electron 
background. 

The difference with the pure magnetic state is particularly striking. We observe a screen- 
ing polarization cloud of s-electrons about the impurity within the range of £ ~ 2 11 ' 6 or r = 
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3.1 x 10 3 (Af/2) .This is about the Kondo length r K . 
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Fig. 6: The net integrated density f^p (x) dx within a distance £ = 2' 

from the d-spin up component of the impurity. Shown are the spin up, spin 

down components as well as the total density and the polarization. The dj-spin 

of 0.93/2 is screened by 0.46 s-electrons within the range of £ ~ 2 11,6 

(or r « 3 x 10 3 A F /2 ). 



The polarization cloud for the ground state of the Kondo impurity (8) is in principle 
identical with the results for the Friedel- Anderson impurity and is discussed in detail in ref. 
[31]. 

5.3 Friedel Oscillation 

Recently Affleck, Borda and Saleur (ABS) [32] showed that the Friedel oscillations due to a 
Kondo impurity are essentially suppressed within a distance of the order of the Kondo length 
rx- They supported their theory by numerical calculations using NRG. Fig. 7a shows the 
universal behavior of their numerical results for many different interaction strengths. Plotted 
is a function F (The actual amplitude is proportional to [1 — F] £~ D where D is the 

dimension of the system). For F = 1 the Friedel oscillation is canceled while for F = — 1 its 
amplitude is doubled). The author could not resist the temptation to evaluate the Friedel 
oscillations with the FAIR method [33]. Fig. 7b shows the FAIR results of the amplitude 
(1 — F) of the Friedel oscillation for two different interaction strengths. This universal curve 
is shown in Fig. 7a as the full blue curve. It agrees well with the numerical results by ABS. 
(The exact form of F (£/£k) is n °t known explicitly). 
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j=0.15 NRG 




-10 -5 5 

Fig. 7a: Reduction F of the iog 2 g/y 

Friedel oscillation is plotted 

versus £/£ K . The amplitude Fig. 7b: The amplitude [1 - F (£/&;)] of 

is proportional to [1 - F (£/£#-)] . the Friedel oscillation for two different 

Calculated by ABS. Kondo energies (from FAIR). 



6 The FAIR States 

At the heart of the FAIR approach are the FAIR states a) Q and b . Therefore a detailed 
discussion of these states is appropriate. The FAIR states are expressed in terms of the 
Wilson states. The latter are described in appendix A and a basic knowledge is required to 
follow some of the arguments of this paragraph. 

A FAIR state is given as aj = J2v a o c t where the states cj, are Wilson states. Now each 
Wilson state cj, represents all the original band states ip k within the energy cell <L U with 
an energy width A v where A v = (( u +i — GO- The composition of the Wilson states cj, is 
cl=Z u 1 ^ 2 ^2 k (p\- This yields for the FAIR state the composition 

t a o t 

v 

This means that the FAIR state consists of the original s-band states <p\ which have the 
amplitude of a^/y/Z^ in the energy cell € u or the occupation |q;q | 2 jZ v . Now we can express 

Kl 2 _ 1 Z - »,2 _ 2 |a 1 2 

Uln 



ZZ„ 1 Ul Z A, 



where Z is the total number of (p k states in the conduction electron band (for one spin) and 
Z/Z u = 2/ A, 
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Therefore the expression p u = |ctQ | 2 / A u represents (besides the factor Z/2) the compo- 
sition of the FAIR state aj in terms of the original band state ip k in the energy cell C u . 

If one plots Pis cLS cL function of energy then one finds a step function because of the finite 
energy width of the cells <£ u of the Wilson states. If one repeatedly sub-divides the energy cells 
(doubling the number of Wilson states) then a smooth function p (£) emerges. This yields 
the FAIR state i n a quasi-continuous energy band. A rather good approximation of 

p (£) can be obtained by interpolation. 

A comparison of Fig. 5 and Fig.6 for the polarization of the magnetic state ^ ms (T) an d 
the magnetic component ^ms (T) of the singlet state ^ $s shows a remarkable difference in 
the polarization about the impurity although the structure of the two states is identical. This 
is particularly surprising since the coefficients A a ^ = A a ,p in the singlet state are roughly 
l/y/2 of the coefficients A a $ of the magnetic state. However, in the singlet state one has 

a finite coupling between ^ms (T) an d ^ms(I)- This shifts the composition of the FAIR 
states a and fej towards small energies. The difference is that the FAIR states in the two 
states have a very different composition. To demonstrate this difference in the compositions 
Pu = \ a o\ 2 / °f a o ( an d &o) are plotted in Fig. 8 as a function of v for the two different 
states. 




Fig. 8a: The density distribution of the p v = 
cLS cL function of cLS cL function of v. 
Fig. 8b: The density distribution of the p v = 

state ^ss 

clS db function of the cell number v 

It would be more natural to plot function of the energy p ((). But for v close to 

iV/2 the width of the energy cells £. v is less than 10~ 6 and any dependence of p{() on the 
energy cannot be resolved on a linear scale. The probability p v increases close to the Fermi 
energy. 

In Fig. 8a for the magnetic state the compositions of and 6q resemble mirror images. 
The function p v has a maximum at small energies of about 10. 

In Fig. 8b the corresponding plot is shown for the singlet state. The weight p u = |ckq | 2 /A„ 
close to the Fermi level is very different for the singlet state and the magnetic state. One 



a o± | / f° r the magnetic state *$?ms 

i 2 

o.q ± | / Aj, for the singlet ground 
Note the difference in scale. 
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observes in the singlet state a maximum of about 400 and the weights p u in aj and 6q are 
essentially identical and not mirror images. The magnetic component of the singlet state is 
in a subtle way different from the magnetic state. 

It should be emphasized again that the two states aj and contain the whole information 
about the many electron states ^ms or ^ss- When Oq and b\ are known the whole bases 
jajjand an d all the coefficients A a> p or A a ^,A a ^ can be reconstructed. 

One important advantage of the FAIR method is that one can modify the size of the 
energy cells <L V after one has completed the numerical calculation of the FAIR states and the 
ground state. It requires relatively little effort to sub-divide the Wilson states. This is an 
important advantage over NRG which can't change the energy width of the Wilson energy 
cells <L V . This procedure is discussed in appendix C. 



A Wilson's s-electron basis 

Wilson [15] in his Kondo paper considered an s-band with a constant density of states and 
the Fermi energy in the center of the band. By measuring the energy from the Fermi level 
and dividing all energies by the Fermi energy Wilson obtained a band ranging from —1 to 
+ 1. To treat the electrons close to the Fermi level at ( = as accurately as possible he 
divided the energy interval (—1 : 0) at energies of — 1/A, — 1/A 2 , — 1/A 3 , .. i.e. ( u = —1/A U . 
This yields energy cells £„ with the range { — l/A^ 1 : —1/A U } and the width A u = C/ — C-i 
= 1/A U . Generally the value A = 2 is chosen. (There are equivalent intervals for positive 
^-values where v is replaced by (N — v) but we discuss here only the negative energies). The 
new Wilson states c* are a superposition of all states in the energy interval (Ci/-i,C>) an d 

have an (averaged) energy (£, + C„_i) /2 = ^, i.e. -|, -|, - A ? .. ? 

(I count the energy cells and the Wilson states from v — 1 to N). This spectrum continues 
symmetrically for positive energies. The essential advantage of the Wilson basis is that it 
has an arbitrarily fine energy spacing at the Fermi energy. 

Wilson rearranged the original quasi-continuous electron states <p\ in such a way that 
only one state within each cell <L V had a finite interaction with the impurity. Assuming that 
the interaction of the original electron states ip\. with the impurity is /^-independent this 
interacting state in <L V had the form 



4 = ZzA/Vz. 

where Z v is the total number of states ip\. in the cell C„ [Z v = Z (( v — Cv-i) /2, Z is the total 
number of states in the band). There are (Z u — 1) additional linear combinations of the 
states (p\ in the cell <L V but they have zero interaction with the impurity and were ignored 
by Wilson, as they are within this work. 
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B Construction of the Basis CLq 5 Qj^ 

For the construction of the state a\ and the rest of basis a\ one starts with the s-band elec- 
trons {cj,} which consist of N states (for example Wilson's states). The c^-state is ignored 
for the moment. 



• In step (1) one forms a normalized state a' out of the s-states with: 



A? 



(10) 



The coefficients can be arbitrary at first. One reasonable choice is «q = 1/yN 



• In step (2) (N — 1) new basis states a\ (1 < i < N — 1) are formed which are normal- 
ized and orthogonal to each other and to aj. 

• In step (3) the s-band Hamiltonian H is constructed in this new basis. One puts the 
state at the top so that its matrix elements are H 0i and H i0 . 

• In step (4) the (N — l)-sub Hamiltonian which does not contain the state Oq is diag- 



onalized. This transforms the rest of the basis 



into a new basis 



a ,a\ \ (but 



keeps the state % unchanged). The resulting Hamilton matrix for the s-band then has 
the form 

/ E(0) Vfr(l) V fr (2) ... V fr (N-l) \ 



Hn 



Vfr(l) 

V>,(2) 



E(l) 








E{2) 









11} 



\ V fr (N-l) ... E(N-l) ) 

The creation operators of the new basis are given by the set |oq, a|| , (0 < % < N — 1). 

The a\ can be expressed in terms of the s-states; a\ = X^Li a i °l- The state uniquely 
determines the other states a\. The state a\ is coupled through the matrix elements 
Vf r (i) to the states a], which makes the state aj an artificial Friedel resonance. The 
matrix elements E (i) and V/ r (i) are given as 

E(i) = J2< £ » a i 

V 
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• In the final step (5) the state % is rotated in the iV-dimensional Hilbert space. In each 
cycle the state is rotated in the ^a^aj^ plane by an angle 9 io for 1 < i < N — 1. 
Each rotation by 9 io yields a new 

ao f = <4 cos 9 io + a\ o sin 9 io 

The rotation leaves the whole basis |aj,a]| orthonormal. Step (4), the diagonaliza- 

tion of the (N — l)-sub Hamiltonian, is now much quicker because the (N — l)-sub- 
Hamiltonian is already diagonal with the exception of the io- row and the io-column 
. For each rotation plane ^aJ,aj o ^ ^ ne optimal a\ with the lowest energy expectation 
value is determined. This cycle is repeated until one reaches the absolute minimum of 
the energy expectation value. In the example of the Friedel resonance Hamiltonian this 
energy agrees numerically with an accuracy of 10~ 15 with the exact ground-state en- 
ergy of the Friedel Hamiltonian [35]. For the Kondo impurity the procedure is stopped 
when the expectation value changes by less than 10~ 10 during a full cycle. 



C Changing the Wilson Basis 

In NRG one usually constructs the Wilson states with A = 2. This means that one uses 
energy cells whose width reduced by a factor of two. NRG is in principle exact for A « 1 
(together with the requirement that one can diagonalize matrices of gigantic sizes). In the 
FAIR method we also begin the calculation with A = 2. When the FAIR state a\ is obtained 

for A = 2 in the basis {cj,} then it is also approximately known in the original basis jv 9 ! j 

(with 10 23 states) 

«o = Ef=i«o4 = Ef=i«oE £l ^i/ \fz~v 

Now one can choose a smaller A, for example A = \[2 and interpolate with good accuracy 
the FAIR state for the smaller value of A [33] . The optimization of the resulting FAIR state 
requires now a relatively short additional numerical iteration. For the calculation of the 
Friedel oscillation a value of A = \[2 = 1. 19 was used. 



D Geometrical derivation of the Friedel ground state 

If the conduction electrons are described by a basis of N states then together with the d-state 
this yields an (N + l)-dimensional Hilbert space $)n+i- The Friedel Hamiltonian is a single 
particle Hamiltonian and possesses in our case (N + 1) orthonormal eigenstates 6] which are 
compositions of the states <4 and the one d state eft. The n-electron ground state is then 
the product of the n creation operator fet with the lowest energy (applied to the vacuum state 
$o)- These n states define the n-dimensional occupied sub-Hilbert space S) n . The remaining 
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(N + 1 — n) eigenstates form the complementary unoccupied sub-Hilbert space S)N+i-n- 111 
the following we treat the creation operators as unit vectors within the Hilbert space. 

Now the vector d of the d state lies partially in the occupied and the unoccupied sub- 
Hilbert space. It has a projection d[ in the occupied sub-Hilbert space S) n and a projection 
d' 2 in the unoccupied sub-Hilbert space S^N+i-n (so that d = d[ + d 2 ). The lengths of 
the vectors d[ and d' 2 are less than one. So we normalize them to di and d 2 with |dj| = I. 
These two vectors are orthogonal (they lie in different sub-Hilbert spaces) and form therefore 
a two-dimensional space. The vector d lies within this plane because 

d = d' x + d' 2 = adi + f3d 2 

The vector perpendicular to d in this plane is the FAIR state a with the composition 

a = /3di — ad 2 

Then the vector di has the form 

di = /3sl + ad 

The vector di can be used as a basis vector of the (N + 1) Hilbert space S)n+i- It lies 
completely within the occupied sub-Hilbert space. Now we divide the occupied sub-Hilbert 
space $) n into the one-dimensional space di and an (n — l)-dimensional subspace @ n _i which 
is orthogonal to d x . This subspace & n -i is also orthogonal to the d state vector d and is 
therefore built only of c u vectors. It can be decomposed into (n — 1) orthonormal basis 
vectors a^. 

Returning to the physics, the ground state can be expressed as 

* F = ((3al + a d ]) Y\a\% 
v 7 i=i 

Similarly the sub-Hilbert space ^Ar +1 _ n can be divided into d 2 and a sub-space &n-u 
orthogonal to d 2 which is therefore also orthogonal to d. &N-n can be expressed in terms 
of (N — n) orthonormal basis vectors a^ (which consists only of vectors c u ). 

The creation operators Oj are not yet uniquely determined. That is done by diagonalizing 

the Hamiltonian Hq in 6„_i and ©Ar_ n . This yields the new basis |aj, 1 < i < N — l|. It is 

straight forward to show that the matrix elements (^a\ \H \ a\,^ for a\ G & n -i an d a\, G &n-u 

vanish as well. (We know the matrix elements of Hp between any state in &N-n and any 
state in @„_i vanishes because the two sub-Hilbert spaces are built from a different sub-set 

of eigenstates of Hp. Therefore (a\ \Hp\ a\^j — if a\ e & n -i an d a\, G &N- n - Since <5 n -i 

and &N-n are orthogonal to eft the d component of the Hamiltonian Hp vanishes anyhow 

and the remaining part ( a\ \Hq \ a\, ) — vanishes for all pairs of i and i'.) 
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